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Recently, clustering of inertial particles in turbulence has been thor- 
_J oughly analyzed for statistically homogeneous isotropic flows. Phe- 

nomenologically, spatial homogeneity of particles configurations is bro- 
ken by the advection of a range of eddies determined by the Stokes 
!— 1 relaxation time of the particles which results in a multi-scale distribu- 

tion of local concentrations and voids. Much less is known concern- 
ing anisotropic flows. Here, by addressing direct numerical simula- 
tions (DNS) of a statistically steady particle-laden homogeneous shear 
OS flow, we provide evidence that the mean shear preferentially orients 

particle patterns. By imprinting anisotropy on large scales velocity 
fluctuations, the shear indirectly affects the geometry of the clusters. 
CN Quantitative evaluation is provided by a purposely designed tool, the 

angular distribution function of particle pairs (ADF), which allows to 
address the anisotropy content of particles aggregates on a scale by 
scale basis. The data provide evidence that, depending on the Stokes 
relaxation time of the particles, anisotropic clustering may occur even 
in the range of scales where the carrier phase velocity field is already 
!_h recovering isotropy. The strength of the singularity in the anisotropic 

component of the ADF quantifies the level of fine scale anisotropy, 
which may even reach values of more than 30% direction-dependent 
variation in the probability to find two close-by particles at viscous 
scale separation. 



1 Introduction 

Transport of inertial particles is involved in several fields of science, e.g. 
droplets growth in clouds, (|TJ EJ, planetary formations, (J3|, or plankton 
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accumulation in the ocean, (j3J). As far as technological applications are 
concerned, inertial particle dynamics is crucial for solid or liquid fuelled 
rockets, injection systems of internal combustion engines or for sediments 
accumulation in pipelines, e.g. d5H6j). Inertial particles differ from perfectly 
Lagrangian tracers due to inertia which prevents them from following the 
flow trajectories. The main effect consists of "preferential accumulation", 
see for instance (J7J Ej). In homogeneous isotropic conditions it amounts to 
the small-scale clustering discussed in a number of recent papers, (jHl El IT0|) . 

In presence of inhomogeneity new features emerge leading to the so- 
called turbophoresis as a preferential accumulation near the boundary in 
wall turbulence (|lip . Under appropriate conditions, particles may achieve 
extremely large concentrations at the wall with a substantial reduction of 
mobility. This turbulence-induced transport and the issuing preferential ac- 
cumulation as been addressed in a number of papers dealing with a variety 
of configurations, from boundary layers to planar channels and pipes, at- 
tacked both from the experimental (1121 [T3l ITS]) and the numerical (|5j EJ fT5|) 
side. Though a complete understanding of the phenomenon is still lack- 
ing, the advection of the particles by the coherent motions in the wall-layer 
is certainly essential, as discussed by (fT6|) . In other words the structures 
responsible for particle accumulation at the wall are the same which sus- 
tain turbulence kinetic energy production in the buffer layer, see also {[Sj [6J) . 
Turbophoresis and small scale clustering are different aspects of the same in- 
ertial particle dynamics. Both phenomena are induced by non trivial phase 
relationships due to quasi-coherent vortical structures. The main difference 
is provided by the characteristic scales, associated with the Kolmogorov 
time unit t„ = tj 2 jv in one case-?] being Kolmogorov length and v the 
kinematic viscosity-and with the larger energy producing time scale in the 
other. Inhomogeneity is essential to have spatial segregation. For instance, 
in the kinetic model presented in (Hip , the spatial transport of particle con- 
centration presents, beside a Fick-like gradient type diffusion component, a 
contribution associated with the spatial variations of turbulence intensity. 
However anisotropy is probably a key ingredient of the process, see e.g. the 
preferential direction of the trajectories of particles approaching the wall 
(|16p . The two features are strongly entangled in wall bounded flows. A 
special flow exists however-the homogeneous shear flow in a confined box- 
which retains most of the anisotropic dynamics of wall bounded flows still 
preserving, on average, spatial homogeneity. 

The flow is bound by a computational box of fixed extension and its inte- 
gral scale grows initially to eventually saturate due to confinement. Target 
of the analysis is the statically steady state with time independent ensemble 
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averages. Similar features are found in the the experimental data by (|17p . 
By using an active grid to generate a flow with integral scale close to the 
transversal dimension of the apparatus, the authors were able to achieve 
confinement from the outstart. As a consequence, the integral scale did not 
increase downstream, see also the discussion in (I18p . 

Our flow shares with the wall-layer streamwise vortices and turbulent 
kinetic energy production mechanisms. In the numerical experiment, this 
corresponds to pseudo-cyclic fluctuations associated with the regeneration 
of streamwise vortical structures. 

Velocity fluctuations are strongly anisotropic at the large scales driven 
by production while, for smaller separations, the classical energy transfer 
mechanisms become effective in inducing re-isotropization. This classical 
issue, see e.g. (fTU]l . has been recently revisited by more complete diagnostic 
tools (e.g. SO(3) decomposition of turbulent fluctuations) able to quantify 
on a scale-by-scale basis the amount of anisotropy in the carrier fluid, as 
discussed both experimentally (j20| I2T1) and numerically (j22l) (see (23) for a 
review) . 

Despite anisotropy of the velocity field is now well understood and the 
carrier fluid shows tendency towards isotropy recovery below the shear scale, 
L s = y/e/S 3 with e the average turbulent kinetic energy dissipation rate per 
unit mass and S the average shear rate, the behavior of particle distributions 
is still not fully explored. 

Anisotropic transport of inertial particles has been recently addressed 
by (|24p , who analyzed numerically the initial transient of the homogeneous 
shear flow, i.e. before saturation of the integral scale occurs, with the 
purpose of modeling unconfined conditions. The focus was mainly on the 
comparison of different particle dynamics models. However, by considering 
particle configurations in orthogonal planes, the authors also discussed the 
anisotropy of particle clusters concluding that particles are most concen- 
trated in the streamwise and least concentrated in the cross-stream direc- 
tion. The same flow was dealt with by (I25[ [26] |2T|) to investigate issues such 
as turbulence modulation in the two way coupling regime or heat transfer 
induced by the disperse phase. 

Purpose of the present paper is the quantitative evaluation of the shear 
induced anisotropy in particles clustering. It is now well known that par- 
ticles respond to the fluid velocity fluctuations in a certain range of scales 
which is determined by their Stokes time. The relevant parameter is the 
Stokes number, ratio of particle Stokes time and flow time scale. In order 
to work with well defined conditions, one needs a shear flow whose charac- 
teristic time and length scales are constant in time. The best candidate, is 
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the statistically steady homogeneous shear flow in a confined box we have 
described above. In this flow, below the shear scale, the velocity fluctua- 
tions tend to recover isotropy. The question is then what happens to particle 
clusters. Do they become isotropic in the smallest scales? For given velocity 
field, how is the geometry of the clusters affected by the relaxation time of 
the particles? 

In fact, the main contribution of the present study is the quantitative ev- 
idence that particle distributions do not necessarily reduce their anisotropy 
at small scales, despite the isotropy recovery occurring in the velocity field. 
Rather clusters anisotropy may even grow below the Kolmogorov length 
where the velocity field is smooth and almost isotropic. As a matter of fact, 
inertia manifest itself in a rather peculiar and unexpected way, and leads, 
under certain coupling conditions, to singular particle distributions which 
viscosity cannot regularize (8). After introducing a suitable observable-the 
angular distribution function- its spherical decomposition is used to evalu- 
ate the relative importance of its different components. The scaling expo- 
nents of the respective singularities show that, under appropriate conditions, 
anisotropy is a leading order effect which may easily persist down to van- 
ishing separations. The data offer preliminary evidence of the anisotropic 
geometry of the fractal support of inertial particle distributions under shear 
flows, thus non-trivially extending results recently achieved for isotropic 
transport. 

2 Methodology 

Concerning the carrier fluid, the velocity field v is decomposed into a mean 
flow U = 5x2 ei and a fluctuation u, see figure [T] for notations. Ro- 
gallo's technique is employed to rewrite the Navier-stokes equations for ve- 
locity fluctuations in a deforming coordinate system convected by the mean 
flow according to the transformation of variables £1 = x\ — Stx2', £2 = 
x 2] £3 = £3; r = t, (|28p . The resulting system is numerically integrated 
by a pseudo-spectral method combined with a fourth order Runge-Kutta 
scheme for temporal evolution, see (|29p . 

The two parameters controlling the homogeneous shear flow are the 
Taylor-Reynolds number Re^ = \/h/{ve){u a u a ) and the shear strength S* = 
S(u a u a ) /e. For the simulations discussed below they are Re a — 100 and 
S* ~ 7, corresponding to a ratio of shear to Kolmogorov scale L s /r] ~ 35. 
Navier-Stokes equations are integrated in a 47r x 2tt x 2ir periodic box with a 
resolution of 256 x 256 x 128 Fourier modes corresponding to 384 x 384 x 192 



4 



p 

V; 



collocation points in physical space due to the 3/2 dealiasing rule. The 
Kolmogorov scale is rj = 0.02 which correspond to K max r] = 3.1 ensuring 
sufficient resoltion at small scales. Actually a well resolved velocity field is 
crucial to minimize numerical errors associated with the interpolation of the 
fluid velocity at particle positions which is necessary to advect the particles 
in the present mixed Eulerian-Lagrangian formulation, see e.g. (|30( I3ip . 

The disperse phase consists of diluted particles with mass density p p 
much larger than the carrier fluid pf, assumed small enough to be modeled 
as material points. At this dilution self-interactions and the back reaction 
on the fluid are negligible, leaving the Stokes drag in the relative motion 
with the fluid as the only relevant force on each particle, (32). Accordingly, 
the equations for particles position x P (t) and velocity vf (t) read 

dx i 
dt 

at Tp 

where Vi(x p ,t) is the instantaneous fluid velocity evaluated at x p {t) and 
T p = Ppdp/ '(l&vpf) is the Stokes relaxation time (d p denotes the particle 
diameter). The particle dynamics is controlled by the ratio of r p to a char- 
acteristic flow time scale, typically the Kolmogorov time scale t v = rf /v, 
i.e. the relevant control parameter is the Stokes number St v = Tp/r^. Par- 
ticle velocities are decomposed as vf = Ui[xf.(t)] + uf where uf denotes the 
particle velocity deviation with respect to the local mean flow of the carrier 
fluid. Eqs. Q can be written as 



to be finally rearranged in Rogallo's computational space as (see also (|24p ) 
dr ^ 



d^ 

— u p , - Stu p 2 5h 



(3) 

(IT Tp 

The particle equations are integrated by the same fourth order Runge-Kutta 
scheme used for the Navier-Stokes equations, with fluid velocity at particle 
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positions evaluated by tri-linear interpolation. The accuracy of the interpo- 
lation scheme may be an issue. To assess its effect on the numerical results, 
we have preliminarily run two different simulations at half the resolution and 
smaller Reynolds number of the cases to be discussed in the main body of 
the paper. The two simulations employ two different interpolation schemes, 
namely a linear and a quadratic one. As always in turbulence, compar- 
isons need to be made in terms of the relevant statistical observables. We 
anticipate that here we deal with the different projections of the Angular 
Distribution Function, to be introduced in § |3j In figure [2] the solid lines 
denote results obtained with the linear interpolation scheme while symbols 
correspond to the quadratic Lagrange polynomials. The difference cannot be 
appreciated on the scale of the diagram, and is always below the statistical 
accuracy of the data. 

Starting from an already fully developed fluid velocity field in statisti- 
cally steady conditions, five different populations of N p = 300000 particles 
each, with Stokes numbers St^ = 0.1, 0.5, 1.0, 5.0, 10.0, are initialized with 
random and homogeneous positions and velocities matching that of the lo- 
cal fluid. Samples for particle statistics corresponding to 120 independent 
snapshots separated in time by 2 S 1-1 are collected after an initial tran- 
sient of 50S* -1 . Discarding the initial transient is crucial to have results 
independent from the rather arbitrary initial state used to initialize the par- 
ticles. We observe that, also under this respect, a statisticaly steady flow is 
mandatory to have well definite experimental conditions, especially in cases 
where the response of particle populations with different relaxation times 
are compared. 

3 Results & discussion 

A visual impression of instantaneous particles configurations is provided in 
figure [3j where slices of the domain in selected coordinate planes are dis- 
played for three different Stokes numbers. The typical particle distribution 
exhibits many voids, strongly correlated with high enstrophy, (jTJ [8|, inter- 
twined with thin "stretched" regions where particles concentrate. Clustering 
is specially manifest near St^ = 1, see the mid panel of the figure in com- 
parison with top and bottom ones which show more even spreading. The 
typical void dimension, as caught by the eye, is larger at our largest Stokes 
number, St v = 10. The ballistic limit, where particles follow their trajec- 
tories with no significant influence from the fluid and the expected spatial 
distribution is homogeneous, is apparently still far away. In the opposite 
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extreme case, passive tracers are recovered for vanishing St„, particles fol- 
low the fluid path, and, again, homogeneity is eventually restored. In fact, 
clustering still takes place, though at smaller scales, for the smallest Stokes 
number we have considered, consistently with theoretical arguments (|10| [1]) 
and numerical simulations (JEJ [9]) for homogeneous and isotropic turbulence 
aimed at explaining droplets growth in clouds. 

Despite of the mentioned similarities with isotropic flows, ours manifests 
specific features compelled by large scale anisotropy. The shear-induced 
orientation is apparent from the bottom- left /top- right alignment of particles 
sheets in the shear plane x—y, see the right panels of figure [3j This behavior, 
clearly visible at St^ = 1.0, is still discernible in the other two cases. Since, 
apparently, the effect is strong and persistent, we are interested in putting 
forward suitable tools to evaluate the anisotropy of clustering. As we shall 
see, this is best done by extending in due form a line of analysis proved 
successful for the isotropic case. 

The main statistical tool is the radial distribution function (RDF) of 
particle pairs g{r) which is a function of radial distance r, see e.g. (|33|) 
where the RDF is dealt with for isotropic flows. The RDF, sometimes called 
correlation function, is defined as 

M * ^1, ,4) 
47rr z ar uq 

where no = 0.5N p (N p — 1) /V$ is the density of pairs in the whole volume Vo 
and N r is the number of pairs in a ball B r of radius r. 

The concept is easily extended to anisotropic cases by considering the 
number of pairs d[i r = v r (r,r)dQ contained in a spherical cone of radius 
r, with axis along the direction f and solid angle d£l, see the sketch in the 
right panel of figure [TJ By this definition the number of pairs in the ball 
B r is N r = J n v r dU, hence dN r /dr = f n dv r /dr dO. We define the Angular 
Distribution Function (ADF) as 



1 dv r 1 
dr no 



2( r > f ) = l2^f— i ( 5 ) 



which retains information on the angular dependence of the distribution. 
The RDF is the spherical average of the ADF g{r) = l/(4-7r) J n g(r,r)dO, 
and it is shown in figure [4] for a few particle populations. 

The behavior of the RDF near the origin, g(r) oc r~ a , is related to impor- 
tant geometrical features of the spatial distribution. Specifically, T>i = 3 — a 
is the so-called correlation dimension of the multi-fractal measure associ- 
ated with the particle density, (I34p . A positive a indicates the occurrence 
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of small scale clustering. Its value is inferred from the slope near the origin 
in the log- log plots shown in figure |4j see the scaling behavior apparent in 
the range r/r] £ [.1 : 1]. From the figure, particles with St v ~ 1 exhibit 
maximum accumulation, i.e. the RDF diverges at a faster rate as r is de- 
creased, see also (24). The solid lines superimposed to the present data 
correspond to the scaling laws extracted at matching Stokes number by ([8} 
in isotropic conditions. The agreement between our data and the isotropic 
ones is remarkable, showing that, even under strong shear, certain features 
of the clustering process may be universal. 

Small scale clustering is not shared by heaviest particles. They show 
instead the saturation of the correlation function to a constant value g{r) ~ 
g* > 1 below a critical scale £ c , see e.g. the open squares in figure [4] with 
a ~ 0. This means that the number of pairs below £ c is proportional to 
volume with an effective density n* = g* uq larger than its overall mean 
value no- By inspection of figure [4] for the heaviest particles (squares), 
the saturation occurs at £ c ~ IO77. The interpretation is that eddies with 
a time scale sufficiently smaller then particle Stokes time do not influence 
the clustering process. The saturation scale should then correspond to the 
size i m i n of the smallest eddies able to aggregate particles. As an order 
of magnitude estimate, our data are consistent with the results given by 
(|35p for the inertial range of isotropic turbulence, where the authors find 
£min/v = {St v / Pmin) 3 ^ 2 with (3 m i n ~ 2. Our data match this estimate also 
for the run at St v = 5, where the saturation occurs close to the Kolmogorov 
scale. For lighter particles, clustering keeps on going below the Kolmogorov 
length and its lower limit cannot be interpreted by arguments which, tuned 
by experiment, are nevertheless taken from inertial range theory. 

Concerning the large scale behavior, heavier particles apparently begin 
to show accumulation at larger scales. Each population is uniform with 
g(r) ~ 1 at very large r. This trend is followed by the lightest particles 
down to 5 T) below which they begin to follow the power-law. The RDF 
starts deviating from the uniform distribution-^?") ~ 1- at a scale £ ma x 
which increases monotonically with the Stokes number, see. e.g. the range 
r/r] £ [20 : 80]. This is consistent with intuition, since the particle Stokes 
time progressively matches the eddy-turnover time of larger turbulent eddies. 
The data again agree reasonably well with the results of (|35p , which identify 
the range of time scales relevant to clustering as I3 max < t p /t(£) < (3 m i n 
where t(£) oc e -1 / 3 /? 2 / 3 is the eddy turn over time at scale I. In terms of 
lengths, clustering starts to occur at min (£o, £ ma x) where £q is the integral 
scale and £ m ax/v = {St v / f3 max ) 3/2 {(3 max — 0.1). Using this estimate, our 
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heaviest particles are expected to begin accumulating at the integral scale 
lo/r] ~ 135 while, e.g., the lightest ones at £ m ax/v — 1) which is not too far 
from the value 5 77 we infer from the plot. 

The RDF quantitatively confirms the overall impression gained from the 
visualizations of figure [3} The strong anisotropy apparent in those plots, 
however, needs a description in terms of the more complete ADF shown 
in figure [5] as a contour plot on the unit sphere for two separations r and 
particles with St„ = 1. The ADF has been normalized with its average on 
the unit sphere, g(r), in order to compare the relative anisotropy content 
of different scales. From figure [5j as separation decreases (i.e. moving from 
the right, r = 35rj, to the left panel, r = 4rj), the normalized ADF exhibits 
preferential clustering in directions consistent with those observed in the 
visualizations of figure [3j 

The ADF allows for a systematic evaluation of anisotropy in particle 
clustering. For given separation r, its angular dependence can be resolved 
in terms of spherical harmonics, 

00 j 

g (r,r)=J2 9jm(r)Y jm (r) . (6) 

3=0 m=-j 

In this notation, the classical RDF g{r) is the projection of the ADF on the 
isotropic sector j = 0, namely 

5oo (r) = g(r) = f g(r, f ) r 00 (r) dfl. (7) 
Jn 

This decomposition amounts in projecting our function on orthogonal sub- 
spaces invariant under rotations. Each sub-space is labeled by the index j 
and it is spanned by 2j + 1 base elements Yj m . Growing levels of anisotropy 
are checked by increasing j, consistently with its geometrical meaning as 
number of zero crossings of Yj m . 

Figure [6] shows the normalized projections on sectors j = 2,4,6, i.e. the 
normalized amplitudes gj m /goo, for the particle population with St^ = 1. 
For given j most modes are negligible and the figure reports only those 
with significant signal level. The (2, —2) mode provides the most significant 
contribution to the anisotropic component of g(r,r). The corresponding 
spherical harmonics, Y2-2, roughly selects the intensity of the signal along 
the principal direction of the mean deformation tensor which corresponds 
to maximum straining and is inclined 45° in the mean flow plane. Y2-2 is 
negative in the first and third quadrant and positive in the others, thus ex- 
plaining the negative sign of 52-2 in figure |6j This description well captures 
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the alignment of most thin particle clusters observed in the mid panel of 
figure [3j 

The signal content rapidly decreases with the order of the sector j and 
the ADF is satisfactorily reconstructed-see decomposition ([6|-using only 
the first few sectors with j < 2 as shown in bottom panels of figure [5j We 
conclude that the probability of finding two particles at separation r = 4r] 
along the main straining direction is 30% higher than in the perpendicular 
direction in the mean flow plane. 

As discussed, the ADF provides a quantitative account of the anisotropy 
induced by the fluid velocity field on the disperse phase. As shown below, it 
can be effectively used to parameterize the level of anisotropy through the 
scales in terms of the Stokes number. 

The plot of figure [7] gives the normalized amplitude of the most energetic 
anisotropic mode in absolute value-|<72-2|/5'oo _ fo r ° ur set of Stokes numbers 
St„ = 10,5,1,0.5,0.1 ranging from heavy to light particles. Focusing on 
the heaviest particles, St v = 10,5, the relative amplitude of the strongest 
anisotropic mode first increases towards the small scales to reach a maximum 
at r ~ £ c . Below this scale the anisotropy level decreases, until the very small 
scales become essentially isotropic. Connecting this result with the previous 
discussion concerning the saturation of the RDF goo , we conclude that the 
heaviest particles show a regular concentration at scales smaller than l c 
where the distribution recovers isotropy. 

Particles with smaller Stokes numbers behave in an entirely different 
way. The anisotropy, as measured by the ratio 52-2/500) substantially in- 
creases to saturate at small scales close to Kolmogorov length. It keeps an 
almost constant value below r/. In other words, the clustering process main- 
tains its anisotropic features even below the dissipative scale for sufficiently 
small Stokes number particles. Remind that the overall clustering process 
described by 500 is here characterized by a singular exponent a. The satura- 
tion observed on the ratio 52-2/900 implies that the dominating anisotropic 
contribution inherits the same behavior, 52-2 oc r~ a . 

4 Final comments 

We have provided evidence that large scale shear induces preferential orien- 
tation on the patterns a disperse phase of small inertial particles forms in 
turbulence. The effect is indirect: The shear imprints anisotropy on velocity 
fluctuations which, in turn, arrange particle configurations in directionally 
biased clusters. 
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Recently it became increasingly clear that the multi-scale nature of the 
velocity field is crucial in explaining most features of particles distributions. 
Typically, a range of eddies exists able to break spatial homogeneity of par- 
ticles configurations. The particles segregate and their pattern shows local 
concentrations and voids on length scales correlated with those of the induc- 
ing eddies. The range of scales of the inducing eddies is determined by the 
Stokes time of the particles, and moves from integral down to Kolmogorov 
length with reducing the relaxation time. For light particles, clustering 
reaches down Kolmogorov scale, leading to a singularity in the radial distri- 
bution function. The exponent of the singularity is a measure of its intensity. 
All these features, originally found in experiments and numerical simulations 
of isotropic turbulence, are also present in our anisotropic fields. They are 
thus generic aspects of the clustering process which seem independent of the 
specific geometry of the forcing. 

At a qualitative level, the specific characteristics of clustering under 
anisotropic advection consists of preferential orientation of the particle pat- 
terns. In order to quantify the new scale-dependent features added by the 
non-trivial geometry of the forcing, we have introduced the concept of an- 
gular distribution function. It can be understood as a generalization of the 
previous radial distribution function, to which it reduces by averaging on the 
unit sphere, i.e. by performing the isotropic projection. This quantitative 
tool has led to our most unexpected finding: The advecting field anisotropy, 
known to be confined to the large scales, affects the singular, small scale, 
clustering process. In fact, anisotropy results into a strong directionality of 
the probability to find a couple of particles at viscous scale separation, with 
30% variations on the solid angle easily observed. 

Technically, for Stokes number order unity, the anisotropic component 
of the angular distribution function diverges at small scales with singularity 
exponent comparable to that found in the isotropic projection (RDF). For 
very small Stokes number, we cannot even exclude that the singularity ex- 
ponent of the strongest anisotropic sector may even exceed that of the radial 
distribution function. Conversely, heavy particles appear to preferentially 
concentrate on finite sized patches endowed with a range of multi-scale and 
shear oriented features, with finest scales more or less evenly and isotropi- 
cally distributed. 

The geometrical properties of patterns of inertial particles differ consid- 
erably from that one could naively guess from velocity fluctuations. Recent 
findings, extending somehow a number of previous results on shear induced 
anisotropy (fT9l [36] [371 EH HOI) show that velocity fluctuations manifest 
two neatly distinct ranges, one dominated by the production of turbulent 
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kinetic energy above the shear scale L s , the other, corresponding to the clas- 
sical inertial range of Kolmogorov theory below. In the two ranges, velocity 
fluctuations display different isotropy recovery rates, a smaller one in the 
production range, a larger one in the inertial transfer range between L s and 
the viscous scale 77, (122 | [2"Tjl . Actually, concerning the velocity field, isotropy 
is matter of fact always recovered at dissipative scales, provided the scale 
separation L s /r/ is large enough, i.e. the local Reynolds number (Taylor- 
Reynolds number, or, equivalently for wall bounded flows the wall normal 
distance in inner units, y + = U\Jt w / p/v, with t w the wall shear stress) is 
sufficiently large. This condition is violated in an essential way only very 
close to solid walls in wall-bounded turbulence (small y + ). 

The anisotropy in the particle configurations depends strongly on the 
properties of the advecting velocity field. However, despite of the isotropy 
recovery of the velocity field, isotropy may never be recovered in the small 
scales of the clusters, as it happens at small Stokes number. Actually, in the 
range from L s to 7] clustering of light inertial particles shows a substantial 
increase of directionality. 

From our results at moderate Reynolds number one can conjecture the 
behavior of the clusters at high Reynolds numbers. In principle the intermit- 
tency of the turbulent field may induce a dependence on Reynolds number. 
In fact, as shown in (JSJ) for isotropic flows, the fractal properties of particle 
distributions depend at most weakly on Reynolds number and strongly on 
the Stokes number. For given geometry of the external forcing, i.e. fixing 
the integral and the shear scale, and at given particle Stokes time, the in- 
crease of the Reynolds number corresponds to increasing the Stokes number 
based on Kolmogorov time, St„. Given the weak dependence on Reynolds 
number, this is somehow equivalent to reading figs. [4] and [7] by successively 
moving from light to heavy particles. Along the process, we infer the sat- 
uration at small scales of the radial distribution function (St^ 3> 1, fig. [4| 
and the small scale isotropy of the clusters (fig. [7]). Clustering is confined to 
the intermediate scales where it shows high levels of anisotropy. This means 
that, in the limit of large Reynolds number, any particles population should 
be organized in anisotropic finite-sized patches which are eventually uniform 
and isotropic in their finest scales. We stress once more that an intermedi- 
ate range of scales always exists, however, where a multi-scale aggregation 
process takes place with a substantial directionality of the clusters. 

On the other hand, at finite Reynolds number, sufficiently small particles- 
small relaxation time-will always show small scale clustering, in the sense of 
a singularity in the radial distribution function. In this case anisotropy may 
persist below Kolmogorov scale, as described by the angular distribution 
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function at small separations for the lightest populations in our simulations. 
This finite Reynolds number effect becomes extremely important in the near 
wall region of turbulent wall bounded flows, where the local Reynolds num- 
ber constructed with the distance from the wall is not huge and the velocity 
field is meanwhile strongly anisotropic. 

The issuing anisotropy of the fine scales of particles clusters will then 
have a significant impact on phenomena of collision, aggregation of dusts 
into larger particles, evaporation/condensation rate of droplets in pipe lines 
and a number of other physically and technologically significant contexts. 

A final comment concerns the extension of the present results to the near 
wall region of wall bounded flows. As shown by recent results on the scale- 
by-scale statistics of the velocity field, e.g. energy transfer across scales, 
spectral distribution of turbulent kinetic energy production, intermittency 
and anisotropy, the homogeneous shear flow reproduces the essential fea- 
tures observed in the wall region, despite of significant differences in the 
large scale geometry of the two systems. However, particles in the wall re- 
gion are strongly affected by turbophoresis which is a predominant effect 
associated with inhomogeneity. In wall flows particle segregation is con- 
trolled by the two concurrent processes of small scale anisotropic clustering 
and of accumulation at the wall. Clearly, the focus of the present paper is 
on the former one, leaving the combined analysis of the two effects for future 
investigations. 
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Figure 1: Left: sketch of the shear flow and nomenclature: the mean flow 
U(y), in the x = x\ direction, is a function of y = X2, with z = x%. For a 
linear mean profile, the shear rate S = dU/dy is constant. Right: sketch of 
the spherical cone of amplitude d£l in direction f . 
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Figure 2: DNS of particle laden homogeneous shear flow (Re\ = 60, S* = 7, 
St v = 1), in a 4-7T x 2tt x 2-7T computational box with 192 x 192 x 96 colloca- 
tions points corresponding to K max r] = 3. The disperse phase is computed 
by using two different schemes to interpolate the fluid velocity at particle 
positions, namely linear interpolation (solid line) and quadratic Lagrange 
polynomials (symbols). The statistical observables shown in the plots are 
defined in § [3| Left panel: <?oo( r )> projection of the ADF in the isotropic 
sector. Right panel: |c?2— 2 1 / <?oo 5 normalized most energetic anisotropic com- 
ponent of the ADF. 
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Figure 3: Snapshots of particle positions for increasing Stokes number, from 
top to bottom: St v = 0.1, 1, 10 respectively. Left column thin slice in the 
y — z plane; right column slice in the x — y plane. The slice thickness is of 
the order of a few Kolmogorov scales. 
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Figure 4: Radial distribution function vs separation, for different Stokes 
number. 
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Figure 5: Top panels: Angular Distribution Function (ADF) giving the 
probability per unit solid angle to find a couple of particle at fixed distance 
|r|. Left panel ADF computed at separation 4ry; right panel 35r/. Data for 
St™ = 1. Bottom panels: Estimate of ADF by using only the isotropic sector 
and the j = 2 sector. 
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Figure 6: Projection of ADF on the different anisotropic sectors of spherical 
harmonics normalized by the projection of the isotropic sector (RDF) as a 
function of separation. Data for St v = 1. 
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Figure 7: Ratio between the most energetic anisotropic sector (2, —2) nor- 
malized by isotropic sector as a function of separation for different Stokes 
number. 
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